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We derive a new factorization relation for the semileptonic radiative decay B nlvy in the 
kinematical region of a slow pion jpV | ~ A and an energetic photon S> A, working at leading 
order in A/mi,. In the limit of a soft pion, the nonperturbative matrix element appearing in this 
relation can be computed using chiral perturbation theory. We present a phenomenological study of 
this decay, which may be important for a precise determination of the exclusive nonradiative decay. 
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I. INTRODUCTION 

Experiments at the B factories are reaching an accu- 
racy level that requires a good understanding of radia- 
tive effects in weak B decays. In particular, real photon 
emission distorts the spectra of other decay products, 
therefore directly affecting the distribution of kinematic 
quantities used to select and count signal events 0, Q ■ 
As a consequence, real photon emission indirectly affects 
branching fraction measurements, leading to systematic 
shifts of the order of a few percent. While radiative ef- 
fects are well studied in semileptonic and nonleptonic 
kaon decays, where chiral perturbation theory is applica- 
ble over the entire kinematical range (see Refs. 01111), 
the situation is less satisfactory in B decays. Here the 
large energy release generally prohibits the application 
of chiral perturbation theory methods over large parts of 
the phase space, and one has to resort to different theo- 
retical tools. 

We study here one of the simplest radiative processes, 
the semileptonic radiative B decays B — > ttIvj. The cor- 
responding nonradiative process, the semileptonic decay 
B — > iriv, is well studied in connection with the determi- 
nation of the CKM matrix element iK^fil- The measured 
branching fractions of this mode are [g 

^ TT-£+iy) = (1) 

(1.33 ± 0.18 ± 0.13) X 10-4 (CLEO) 
(1.38 ± 0.10 ± 0.18) X 10--* (BABAR) 
(1.48 ± 0.20 ± 0.16) X 10-^ (BELLE) 

which contribute to the world average B{B^ 
TT-e+iy) = (1.35 ± 0.08 ± 0.08) X 10-4 0. 

The hadronic dynamics of the radiative semileptonic 
decays is complicated by the presence of many energy 
scales. However, in certain kinematical regions a model 
independent treatment is possible. In the limit of a soft 
photon. Low's theorem "s*] fixes the terms of order 0{q^^) 
and 0{q'^) in the expansion of the decay amplitude in 
powers of the photon momentum q. These contributions 
have been computed in 9] in kaon semileptonic radiative 
decays Kg^. Higher order terms in q have been more 
recently obtained using chiral perturbation theory 0, |Q . 
A similar approach is possible in B decays with the help 



of the heavy hadron chiral perturbation theory jlOl [III 
01 , as long as the pion and photon are sufficiently 
soft. This is the case in a limited corner of the phase 
space, corresponding to a large dilepton invariant mass 

The radiative semileptonic decay B — > "fniDi at low 
produces in general either an energetic pion or an 
energetic photon (see Fig. 1). In this pap er we use soft- 
collinear effective theory (SCET) [11 [ll El El -[Hill 

methods to prove a factorization relation for the decay 
amplitude in the kinematical region with an energetic 
photon and a soft pion. This can be written schemat- 
ically as 

A{B ^ irei'j) = H (g) J (g)S{B^TT) + 0{A/Q) (2) 

with A ^ ^QCD and Q ^ {E^,mi,}. Here H and J 
are hard and jet functions, calculable in perturbation 
theory. S{B — > tt) is related to the matrix element of 
a nonlocal soft operator which describes the B ^ tt 
transition. This factorization relation generalizes to any 
process B ^Mtv [M = p,ui, - ■ ■ a light soft meson) 
with the appropriate soft function S{B M) replac- 
ing S{B — > tt). In this work we use heavy hadron chiral 
perturbation theory to compute the soft matrix element 
S{B tt), and express it in terms of one of the B meson 
light-cone wave functions (j)-^{k+) The convolution 

of this function with the jet function J in Eq. |(2Jl is sim- 
ilar to the one appearing in a factorization relation for 
the radiative leptonic decay B ^ ^ii^ [21II2I I2ll2i| . 

The paper is organized as follows: after fixing the no- 
tation in Section m we describe the main results of the 
paper (factorization formula for hadronic form factors) in 
Section Ulll We give a proof of the factorization formula 
in Section llVI In Section we explore the phenomeno- 
logical consequences of our results, computing the ra- 
diative branching fractions with appropriate cuts on the 
photon and pion energy, as well as various differential de- 
cay distributions. The factorization results are compared 
with the ones obtained by extrapolating Low's amplitude 
to the region of an energetic photon. We collect our con- 
clusions in Section IVTI An appendix contains the expres- 
sion of the radiative amplitude in Low's limit of a soft 
photon. 
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II. AMPLITUDE DEFINITIONS 

We consider in this paper the radiative semileptonic B 
decays 



B{p) ^ Tr{p')e {Pe)i^{Pu)l{q,£) 



(3) 



mediated by the combined semileptonic weak Hamilto- 
nian 



Hw — 



AGf 
V2 



(4) 



and the electromagnetic coupling of quarks and charged 
leptons. 

The full amplitude for B{p) — > 'K(p')e~ {pe)v(jpv)^{q, s) 
contains two terms, describing the photon coupling to the 
hadronic system and to the charged lepton, respectively 



A 



eGp 

71" 

F-'{p,p') 



2q-p 

We denoted here 

F^p') - {7rip')\uj,b\Bip)) , 
and the hadronic tensors T^'^ are defined as 



(5) 
(6) 



T/, = ^ j d4^e^.-(^(p')|T , J.(0)} \B{j>)) (7) 

with J™ = eqQlp,Q the quarks' electromagnetic cur- 
rent and: 



T, 



V 



flu 

A 



= K = U7^& (8) 

• Ju= A^^ W7i^75&- (9) 
The hadronic tensors satisfy the Ward identities 

q^T^, = -F,{p,p') (10) 

q'T^. = . (11) 

The most general form for the vector amplitude com- 
patible with the Ward identity can be written in terms of 
4 invariants V1-4. After multiplication with the photon 
polarization, this is written as (W — Pe + Pu) 



Vi{el 



W 



Hp' ■ £ 



'"^ q.W 
{p'-q){e*-W) 



e* ■ W 
' q-W 



q ■ W 
F.ip,p'). 



qu) (12) 
){V2q,. + V^W, + Vip',) 



The axial amplitude can be written in terms of 
another four invariants Ai_4. They can be chosen as 

e*%t = ^e^,p„e*^'{A^p'pq, + A2qpW^) (13) 
+ ienXpae*''p'xqpW„{A:iW^ + A^p'^) 



The invariant form factors Vi_4 and Ai_4 depend 
on 3 independent kinematical variables, which can be 
chosen as W'^ , E-^ , E-^ with Et^ and E^ defined as the 
pion and photon energy in the rest frame of the B me- 
son. For given W^^, the possible final states in the de- 
cay can be represented as points in the Dalitz plot for 
{Ej,E^). We show in Fig. 1 the Dalitz plot for ^ 4 
GeV^. The dilepton invariant mass can take values 
= [0, (tob - m^)2] = [0, 26.5] GeV^. 

The form factors in B (quark content bq) and B (bq 
quark content) are related by charge conjugation. The 
transformations of the currents under charge conjugation 
are 

Cq-f„bC^ = -Irf^q, Cqj^j^bC^ = hi^lBq (14) 

We adopt everywhere in this paper the following phase 
conventions for the meson states 

(tt^ . 7r° , TT^) — (ud, —r={uu — dd), du) (15) 
v2 

{B~,B°) = {bu, bd) , (S°, B+) = (bd, bu) . (16) 

With this convention the states transform as Ctt^ = 
TT- ,CS° = B" ,CB~ = B+. In addition, the tt^ is C- 
even and the photon is C-odd. The semileptonic radiative 
_B — !■ TT form factors transform as follows under charge 
conjugation, independent of the flavor of the spectator 
quark 

Vi-4{B TTje-De) ^ Vi-i{B TT-fi+iyg) (17) 
^1^4(5 TTje^Pi) = -Ai^4{B TT-yg+iyi) (18) 

The vector form factor F,y{p,p') is C-odd. Finally, the 
form factors depend also on the flavor of the light spec- 
tator in the B meson. For simplicity of notation, we chose 
to omit this dependence, but it will be specifled explicitly 
whenever required. 



E^(GeV) 




0.5 



1.5 



E^(GeV) 



FIG. 1: The Dalitz plot for B -f-Klv at = 4 GeV^ Our 
computation of the form factors is applicable in the lower right 
corner of the Dalitz plot, with an energetic photon and a soft 
pion. 
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III. RESULTS 

We describe here the main resuhs of this paper. In the 
kinematical region h and i?^ ~ A (hard photon 

and soft pion), the form factors for B iriv^ are given 
by a factorization theorem. At leading order in A/Q with 
Q {rribjEj}, these factorization relations are the same 
for both B~ ,B^, up to an isospin factor. The predictions 
for the vector form factors in 7r+£~P7 are 



V2 



2euC'f \E^) J dk+J{k+)Si{k+,f,0 (19) 
^{2Ci^\E,) + ctHE,) + 2C!,''\E,)) 



X / dk+J{k+)S2{k+,t^,0 



(20) 



2eu ^(v) . 



n-W 



C),'{E^) / dk+J{k+)S2{k+,t\0 (21) 



(22) 



and the corresponding relations for the axial form factors 
are 



A4M2 



^^{2&^\e,) + &^\e,) + 2&^\e,) 

X j dk+J{k+)S2{k+,t\0 



(23) 



q-W 



c'i^'He^) J dk+J[k+)Si{k+,t\0 (24) 
&^\e,) (25) 



q-Wn-W 

X j dk+J{k+)S2{k+,t\0 



(26) 



the photon momentum. 

The hard functions C\"}^{E^) are Wilson coefficients in 
the SCET, and are calculable in perturbation theory in 
an expansion in as{Q) as C^'^ = 1 + 0{as{mi,)) , c'^l^ = 
0{as{rnb)). Their one-loop expressions can be found in 
Ref. [13 ■ The jet function J{k^) is calculable in an ex- 
pansion of as(\/ QA) and its one-loop expression is given 
below in Eq. (|36|l . Finally, S'i_2(fc+, i^, C) ^re nonpertur- 
bative soft matrix elements, defined below in Eqs. (|38|l . 
(123. 

The results (|19|l - (|26|l show that, at leading order in 
A/Q, the eight form factors parameterizing B — > ^i:ev 
are given by only two nonperturbative parameters and 
thus satisfy symmetry relations. 



The soft functions depend in general on k+, t'^ = 
{p — p')'^, ( = n ■ p'/n ■ p, and are calculable at leading 
order in chiral perturbation theory. In the soft pion limit 
they can be expressed in terms of one of the B meson 
light-cone wave function and are given by 

S.ik,,t^O = -^^?(Mfl+5^^)(27) 



4/. 



vp' + A 



1 



vp' + A 



(28) 



where v — p/ruB, A — tjib' — ruB — 50 MeV, 63 — n~v, 
and g is the axial effective coupling defined in Eq. H45|l. 
To this order, the dependence on nonperturbative dy- 
namics of the B meson enters only through one quantity, 
the convolution of the light-cone B wave function 0:^(fc+) 
with the jet function J(fc+): 



Vi, cx / dk+J{k+)cl)+{k+,fi) . 



After multiplication by the appropriate tensor structures, 
V4 and A4 give subleading contributions in A/Q to the 
hadronic tensor. The corresponding form factors for 
B~ ^Tr°£~D are multiplied by a factor 1/V2. In these 
formulas — q^/E^ is a light-cone unit vector along with n'^ 



(29) 



A similar convolution appears in the factorization rela- 
tion for the radiative leptonic B decay B —^ jeD, the only 
difference being that in this case the argument of the jet 
function is time-like {+2E~fk+) as opposed to spacelike 
in the radiative leptonic case. As a consequence the form 
factors Vi and Ai are in general complex, receiving strong 
phases of ~ as{^/2E^) via the loop corrections to the 
jet function J. 



IV. FACTORIZATION PROOF 

The proof of the relations Eqs. 119|) - (|26ll makes use of 
the soft coUinear effective theory (SCET) J^EllllIll 
ITM Hgf . The basic idea will be to match the correlator of 
currents in Eq. Q onto soft SCET operators, and then 
in a second step to construct the chiral representation of 
the resulting soft operators. The i? — > tt matrix element 
of the correlator in Eq. (|7|) will be computed in the low 
energy chiral perturbation theory. 

We start by setting up the kinematics of the pro- 
cess. We introduce two light cone unit vectors n, n 



,n ■ n 



2. The dilepton momen- 



tum W is chosen along n with light cone coordinates 
W = {n-W,n-W, W±) (such that n ■ W ^ n ■ W,W±). 
The photon is emitted along n^, such that its light-cone 
coordinates are q = {0,n ■ q,0), and the pion momentum 
is soft, with components Ptt ^ A. 

We are interested in the correlator of the weak semilep- 
tonic current with the electromagnetic current 

= * / d'xe^^--T{iq^Mm,f^'^ix)} (30) 



The matrix elements of this correlator define the hadronic 
tensors in the decays of interest Eqs. O. The flavor de- 
pendence of the form factors follows from the transforma- 
tion properties of this correlator under flavor SU(3). This 
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is given by3(8)8 = 3 + 6 + 15, so that the _B — > tt transi- 
tions are in general described by three unknown reduced 
matrix elements. However, at leading order in A/E^, 
only one amplitude contributes. To see this, consider 
the two types of contributions to the correlator T^i^(g), 
corresponding to the photon attaching either to the cur- 
rent quarks (u and b), or the spectator quark. At lead- 
ing order in A/Ej, only the first amplitude contributes, 
while the spectator effect is power suppressed (see be- 
low). The dominant amplitudes transform as a 3 under 
flavor SU(3), so they are independent of the flavor of the 
spectator quark. 

The semileptonic 6 — > u current is first matched onto 
SCET/, an effective theory containing ultrasoft quarks 
and gluons with momenta ps A, and coUinear modes 
with virtuality ~ QA moving along ~ n^^ . Ne- 
glecting light quark masses, the V — A current is matched 
at leading order in A = A/Q only onto structures contain- 
ing one left-handed light quark field |15(| 

q-/^PLb = Cl") {Lu)qn,ul^PLby (31) 

+ {ci''\Lo)v^ + Gl'\io)-^)q^.Mv + 0(A) 
n ■ V 

where c|^''(a') are Wilson coefficients calculable in per- 
turbation theory as an expansion in as{Q). 

We match first the correlator Eq. ^ onto SCETj 
operators. The weak current is matched as shown above 
in Eq. (|3f |l . The light quark e.m. current coupling g^„^„ 
of an energetic (n— coUinear) photon to a ultrasoft quark 
is more complicated, and it contains both a local and a 
nonlocal term 

= e, r + r{ J« , J d^xiC^^^ [x]} (32) 

where J^^ describes the coupling of a rt— collinear pho- 
ton to 71— collinear quarks. This scales like 0(A^^) and 
it appears in time-ordered products with a subleading 
ultrasoft-collinear Lagrangian. It is given by 

J« = e,q,,^,^*^[wUl?^^^^^^ (33) 

Finally, we match the SCET/ expression of the cor- 
relator H30|l onto SCET//. After soft-coUinear factoriza- 
tion the coUinear quarks and gluons with virtuality 
~ E^A decouple from the soft modes. This leads to 
the final form of the time-ordered product in Eq. I|3U|) in 
the effective theory 

e*^'T^M ^ -C^i^ e„ j dk+J(k+)Oi,{k+) (34) 
-[C^^^v, + (01"^+ Ci"^)^]e„/dfc+j(fc+)02(fc+) , 

Ti • V J 



with the two soft operators given by nonlocal quark bi- 
linears 

Oi.(fc+) = /^e-i^'=^g-(A^)r„(A|,0)^l|7,^P/.6„(0) 

02(fc+) = /^e-i^'^^g(A|)r„(A|,0)^l|p/^6.(0). 

(35) 

Here Yn{x,y) = Pexp{ig dan ■ A{an)) is a Wilson 
line along n^^. 

The jet function J describes the effects of the coUinear 
quarks and gluons with virtualities p^ ~ E^A. It can 
be written as a vacuum expectation value of SCET/ 
coUinear fields [2^, and as such it is calculable in per- 
turbation theory in an expansion in as{\/QA). The 
explicit result to subleading order in as is (with L — 
log[{~2E^k+ - ie)li?)) 

Finally, the nonperturbative physics is contained in 
the matrix elements of the soft operators Oip(fc+) and 
02{k+) in Eq. 

An alternative derivation is possible, wherein one 
matches directly from QCD onto SCET//. In this ap- 
proach I^jE^Ii the jet function J together with the hard 
coefficients CI are obtained as the Wilson coefficient 
in the QCD — > SCETn matching. We will use here the 
2-step approach, which has the advantage of allowing a 
transparent connection between the SCET/ Wilson coef- 
ficients in H31|l and the Wilson coefficient of the SCET// 
operator. 

The discussion above parallels closely the proof of fac- 
torization in i? ^ ^iv decays in Ref. 23] , except that we 
were careful to keep it independent on the final hadronic 
state, which is only assumed to be soft. For example, we 
included also the operators in the second term of Eq. 1)34(1 
which do not contribute to leptonic radiative B decays. 

To complete the proof of factorization one has to show 
that no other operators can appear at leading order in 
Eq. 134|) . Since soft-coUinear factorization in SCETi pre- 
serves the form of the soft operators, it is clear that the 
T-products of the operators contributing to the weak cur- 
rent Eq. H31|l . and to the e.m. current Eq. I|32|) . will 
always reduce to combinations of Oi^ and 02- Photon 
couplings to the spectator quarks require at least 2 in- 
sertions of the soft-coUinear Lagrangian , so they are 
power suppressed. Finally, no operators containing ad- 
ditional soft quarks and gluons appear at leading order; 
the proof is similar to that for B ")ev, to which such 
operators would also contribute. For a recent study of 
these subleadingcontributions in the radiative leptonic 
decay, see Ref. |2q . 

Taking the B — s- vacuum matrix element of the 
SCET// operators in Eq. H34|) gives directly the hadronic 
tensor relevant for the radiative leptonic B —f ^iv. This 
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can be evaluated explicitly using the B —f matrix ele- 
ment of the soft operator Oi^{k^). The latter is parame- 
terized in terms of the B meson light-cone wave functions, 
defined as H 



dX 



]q.{\-)Y^{\,QnmB{v)) ^ 



(37) 



i f 1 + 



With this definition, the wave functions are normalized 
as J dk^(j)±[k^) = 1. Using the soft matrix element 
Eq. H37|l into Eq. 1)34(1 reproduces the leading order fac- 
torization relation for the radiative leptonic decay form 
factors 

Proceeding in a similar way, by taking the matrix el- 
ement of the operator identity Eq. (|34|l for the B 
TT transition (with tt a soft pion) gives a factoriza- 
tion relation for the semileptonic radiative decay B — )■ 
TTjev. This requires the introduction of the soft func- 
tions S'i,2(fc+, t^, C) defined by (with t = p — p' and 
C = {n-p')/{n ■ v)) 

{iT{p')\0,,{k+)\B\v)) (38) 

= -{9t. + ^^t^)CSi{k+,eX) 
{^{p')\0,[k+)\B\v)) (39) 

= -[gtu + ^4uVuP'tj.S2{k+,t' ,0 . 

Similar nonperturbative functions appear in the theory of 
off-diagonal hard-scattering exclusive processes , and 
are known as GPDs (generalized parton distributions). 
See Refs. '231 for recent reviews. In the context of 6 ^ u 
decays, they were considered previously in Ref. |29j, and 

Si is directly related to the function JF^^-* (x, g^) intro- 
duced in that paper. 

The soft functions Si ^2 satisfy several model- 
independent constraints. T invariance and complex con- 
jugation constrains them to be real. Their zeroth mo- 
ments over fc+ are related to the usual J5 — > tt form fac- 
tors (see the Appendix for the definitions of the form 
factors) 



1 1 

dk+Si{k+, t\ = --n ■ pUit') + -f2{t^)n ■ t 



1 

Ak+S2{k+,t\0 = -n-pfrit^). 



(40) 



while higher moments are related to the corresponding 
form factors of higher-dimensional operators 



dfc+(fc+)^5i(A:+,i^C) 
= -{n{p')\q{-n-iDf'^PLK\B{p)) 



pI^ / dfc+(fc+)^52(fc+,^^C) 
Jo 



(41) 



(42) 



{nip'M-n■^Df'L^b^\Bip)) 



B 



FIG. 2; Chiral perturbation theory graphs contributing to 
the B TT matrix element of the soft operators Oifi,{k+) and 
02(fc+) (represented by filled circles) appearing in the factor- 
ization relation for B ~* yK£i>. The photon is not shown. 



Inserting the definitions of the soft functions Eqs. I|38|l . 
(|39|l into the operator expression for the hadronic tensor 
Tf^u Eq. H34|l . one finds the leading order factorization 
relations hsted in Eqs. (fll^ - f^ . 

Similar relations can be given for other semileptonic ra- 
diative B decays, such as B — > pjeP and B — > K^e~^e~ , 
in terms of appropriately defined soft functions. In or- 
der to make detailed predictions, some information is re- 
quired about the soft functions Si ^2- In the next section 
we show that in the limit of a soft pion, chiral symmetry 
fixes the functions 5*1^2(^-1-; ti C) in terms of the B meson 
light-cone wave function defined in Eq. (|37|) . 



A. Soft pion matrix elements of Oifi,02 

The i? — !■ TT matrix elements of the soft operators 
Oi^t(fc+) and 02{k+) are constrained by chiral symmetry 
and heavy quark symmetry. The appropriate theory tool 
to derive these constraints is heavy hadron chiral per- 
turbation theory [Ulil El [3. This effective theory 
describes the interaction of Goldstone bosons with them- 
selves and with heavy hadrons. The Goldstone bosons 
are contained in the nonlinear field ^ ~ e^^/f (S = ^2-^^ 
where / is the pion decay constant / ~ 131 MeV, and M 
is the usual SU(3) matrix 



M = 



TT 

K- 



->° + 76^ X" I (43) 
^0 2 



V6 



?7 



The velocity-dependent heavy meson fields {Bv,B*) 
are grouped into the supcrfield 



(44) 



where a is a SU(3) flavor index labelling the flavor of the 
spectator quark in the B'-*\ 

The interactions of these fields are described by the 
effective Lagrangian |ii|ll|ll|ll 

C = ^Tr (9^Ec)^St) _ ^Xr H'^v^d^Ha 

+ |Tr H^H,^,^, [Cta--^ - ^a-^t] + (45) 
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The symmetries of the theory constrain also the form of 
operators such as currents. For example, the left handed 
current L'^ = qa'y^PhQ in QCD can be written in the low 
energy chiral theory as [lol llll IT^ 



(46) 



where the ellipsis denote higher dimension operators in 
the chiral and heavy quark expansions. The parameter 
a is obtained by taking the vacuum to B matrix element 
of the current, which gives a = fB^/rns (we use a non- 
relativistic normalization for the states as in [isjl. 

The chiral realization of the nonlocal soft operators 
appearing in l|34|) is constructed in a similar way |20| . 
Both operators Oi_2(fc+) contain only the left chiral light 
quark field, so they transform as (3l, 1_r) under the chiral 
group. Making the flavor index explicit, they are both of 
the form (see Eq. 1(23)) 



(47) 



In analogy with the local current (|46ll we write for 
Op(fc_(.) in the chiral theory 



OUk+) = -TT[aL{k+)PRrH4j, (48) 

where the most general form for aL^k^) depends on four 
unknown functions ai{k+) 



elements to leading order in the chiral expansion: 

{7T{p')\0i^{k+)\B{p)) = ^fsmBMk+) (53) 



63 -p 



{^{p')\02{k+)\B{p)) = -^fBmBct>+{k+) (54) 



^ V ■ p' + l\ 

where the unit vector 63 is defined as 63 = ^( n — n). Con- 
fronting these expressions with Eqs. (|38|I - H39|I one obtains 
the results for 5i.2 reported in Eqs. (|?7|l - lf^ . 



V. PHENOMENOLOGY AND NUMERICAL 
STUDY 

In this section we present some phenomenological ap- 
plications of our factorization formula, with a twofold 
intent: (i) give predictions for quantities such as photon 
spectrum and the integrated rate with appropriate cuts; 
(ii) compare our effective-theory based results with sim- 
plified approaches to i? — > 7rej>7, which are usually im- 
plemented in the Monte Carlo (MC) simulations used in 
the experimental analysis. All numerical results reported 
here refer to the neutral B decay 7r"'"e~i?e7- 



aL{k+) = ai + a27^ + a-s-p + ^a4['^,-fi\ (49) 



Factorization results 



The heavy quark symmetry constraint H-p = —H reduces 
the number of these functions to two. Taking the vacuum 
to B meson matrix element fixes the remaining functions 
as 



aL{k^ 



-fBV^W+{k+)+i<l)-{k+)] (50) 



where (/)±(fc+) are the usual light-cone wave functions of 
a B meson, defined in Eq. H37|) . 

We are now in a position to compute the i? — s- tt matrix 
element of the soft operators in Eq. There are two 

chiral perturbation theory contributions shown in Fig. 2. 
Considering the example of the operator Oi^, the two 
diagrams are given by 



We present here results for the photon energy spec- 
trum, the distribution in the electron-photon angle, and 
integrated rate with various kinematical cuts, forced 
upon us by the theoretical tools used: 

1. An upper cut on the pion energy in the B rest frame 

< E'^^^. This restriction is required by the ap- 
plicability of the chiral perturbation theory compu- 
tation of the soft matrix elements of Oi^.2. 

2. A lower cut on the charged lepton-photon angle in 
the rest frame of the B meson, 9^^ > 6'™". This is 
required in order to regulate the coUinear singular- 
ity introduced by photon radiation off the charged 
Icpton leg in the massless lepton limit. 



la = —fBmB<p+{k+){gj^^+iel^^)e*^ 



h = 'i^.fBmB(l)+{k+){g^^ + ief^^)e*^ 



-i{na - Van -v) . 2^ . 



(51) 
(52) 



Adding them together and performing a similar calcula- 
tion for 02(^+) gives the final result for i? — > tt matrix 



3. Finally, to ensure the validity of the factorization 
theorem, we require the photon to be sufficiently 
energetic (in the B rest frame) E^™ > Age d ■ We 
1 GeV. 



choose E^"^ 



In our evaluation we use the transition matrix ele- 
ment in Eq. |SJl, with the factorization expressions for 
the hadronic amplitudes T^^. Working to leading order 
in Qfs, the only surviving form factors are Vi^2 and Ai,2 
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and the nonperturbative input needed is the first inverse 
moment of the light-cone B wave functions <j)^{k^): 



(55) 



Moreover, the tt^£ v transition matrix element 

FviP^p') is parameterized as 

F,{p,p')^f+{E^)ip + p'),, + f^{E^)ip-p'), (56) 

We adopt for this matrix element the leading order 
HHxPT result [13, [H El [H 



fsmB 



1 



2U E^ + A 



(57) 



The corresponding result for B~ —> t:^1~D is multiplied 
by a factor 1/ \/2. We summarize in Table^the numerical 
ranges used for the various input parameters [s^, HI] . 

The phase space integration with the cuts described 
above was performed using the Monte Carlo event gen- 
erator RAMBO l^^l- The factorization results for the 
photon energy spectrum and the distribution in cos Oe-y 
are shown in Fig. |3| (solid lines), using the central and 
extreme values for the input parameters as reported in 
Table m The shaded bands reflect the uncertainty in our 
prediction, which is due mostly to the variation of A^. 
For = 1 GeV, E'^'^'' = 0.5 GeV, and 0™" = 5°, the 
integrated radiative branching fraction predicted by our 
factorization formula is (up to small perturbative correc- 
tions in Q!s(mb) and as{^/rniK)) 

..e.7(fact) = (l.2 ± 0.2(g) 1§J(Ab)) x lO"*^ 

\Vub\\' f fB 



0.004 



200 McV 



(58) 



Apart from the overall quadratic dependence on Vub and 
/b, the main uncertainty of this prediction comes from 
the poorly known nonperturbative parameter As. 

The radiative corrections to these predictions are 
known to one-loop order. Their numerical impact has 
been computed in the radiative leptonic decay case, 
including a resummation of large logarithmic terms 
log^{2E^/li) El m m m and log^{2E^k+/ii'^) 
Their effect can be significant, and was found to reduce 
the tree level result by about 30%. We do not include 
them in our results, in view of the large uncertainty from 
the input parameters, and leave a study of their effects 
for future work. 

Finally, let us mention that radiative rate, photon 
spectrum and angular distribution could be used in the 
future to obtain experimental constraints on A^, given 
their strong sensitivity to it. This requires, however, a re- 
liable estimate of the impact of chiral corrections, 0{as) 
corrections, and the uncertainty in other input parame- 
ters. 



TABLE I: Input parameters used in the numerical computa- 
tion f3^,f3llpl. 



fB 


(200 ± 30) MeV 


u 


131 MeV 


As 


(350 ± 150) MeV 


9 


0.5 ±0.1 


IKbl 


0.004 







TABLE IL Branching fractions of the cut B ■nev'^ decay 
corresponding to the factorization amplitude and the TBI, IB2 
cases defined in the text, for central values of the input pa- 
rameters reported in Table 



cuts 


lO'' Br(fact) 


10^ Br(IBl) 


lO'' Br(IB2) 


E., > 1 GeV 








< 0.5 GeV 


1.2 


2.4 


2.8 


Oej > 5° 









B. Comparison with simplified approaches 

In the Monte Carlo simulations used in most experi- 
mental analyses of B decays, simplified versions of the 
radiative hadronic matrix elements are implemented. 
These rely on universal properties of radiative amplitudes 
following from gauge invariance, which allow one to de- 
rive in a model- independent way the terms of order 
and of the radiative amplitude (see Appendix fXl for a 
summary of these results, based on Low's theorem 8]). 
These simplified amplitudes do not take into account 
hadronic strucure effects and their validity is restricted 
in general to emission of non-energetic photons (whose 
wavelength does not resolve the structure of hadrons in- 
volved in the process, E^ <C A). They are nonetheless 
used over the entire kinematical range, for lack of bet- 
ter results. Our interest here is to evaluate how these 
methods fare in the region of soft pion and hard photon 
(certainly beyond their regime of validity), where we can 
use the QCD-based factorization results as benchmark. 

We use two simplified versions of the B neDj ma- 
trix element, which we denote IBl and IB2 (Internal 
Bremsstrahlung 1 and 2). They arise as special cases 
of Low's amplitude that we report in Appendix 1X1 

IBl: This is the 0{q^^) component of Low's amplitude 
(first two term in the first line of Eq. (jAip ). given by 
the product of the non-radiative amplitude with an 
appropriate radiative term. The PHOTOS Monte 
Carlo |34l] implements this form of the radiative 
amplitude. 

IB2: This corresponds to using Low's amplitude of 
Eq. (|A1|) and setting to zero the derivatives of non- 
radiative form factors (i.e. just keep the first line 
in Eq. HAl|l ). This type of amplitude has been used 
in a series of papers on Ki^^ by Ginsberg [s^l- It 
is often used in B-physics applications by simply 
rescaling the appropriate meson masses. 
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E., (GeV) 

FIG. 3: Left panel: DifFerential rate dV/dE-y, as a function of the photon energy in the B-meson rest frame. The normalization 
is given by F = 10^ F~o x F(i3 — > ireu^). Right panel: Differential radiative rate dT/d{cos9ej), as a function of the cosine 
of electron-photon angle in the B-meson rest frame. The shaded bands and the superimposed continuous lines represent the 
uncertainty in the factorization result (FACT) induced by the parameter As, assumed to range between 200 (upper line) 
and 500 (lower line) MeV, with all other parameters equal to the central values of Table HI The dark(light)-grey dashed line 
corresponds to the IBl (IB2) results (see text for their definition). 



In both cases the only input needed to perform the com- 
putation is /+(i?^) (sec Eq. (|57|l ). 

The results for the photon energy spectrum and the 
distribution in cosOej are shown in Figs. |3| (IBl: dark- 
grey dashed lines, IB2: light-grey dashed lines). The 
integrated branching fractions are listed in Table (to- 
gether with the factorization prediction) for central val- 
ues of the input parameters given in Table HJ The IBl 
and IB2 branching fractions scale very simply with the 
input parameters, being proportional to {\Vub\ g fs)'^ ■ 

As can be seen from the comparison of differential dis- 
tributions and BRs, the two cases IBl and IB2 lead to al- 
most identical predictions, not surprisingly (and in agree- 
ment with previous claims that PHOTOS-type distribu- 
tions are in agreements with Ginsberg's ones P, l33|). 
Both IBl and IB2 amplitudes, however, produce results 
quite different from the factorization ones, as can be seen 
from Fig. O and Table ITTl For central values of the input 
parameters reported in Table the IBl and IB2 predict 
a radiative BR roughly a factor of two larger than the 
factorization result. Moreover, the shape of the distribu- 
tions in and cos9ej are quite different, as can be seen 
from Fig. |31 

There is only one corner of parameter space where the 
integrated rate from factorization gets close to the IBl- 
IB2 result, namely for As 200 MeV. Even in this case, 
however, the distributions in photon energy and electron- 
photon angle clearly distinguish the two scenarios. This 
is explicitly shown in Fig.|21 where the uncertainty in the 
factorization result induced by As is explicitly displayed 
through the shaded region. 

We conclude that, in the kinematical region of soft pi- 
ous and hard photons, the QCD-based factorization pre- 
dictions for the radiative decay are quite different from 
the ones based on the simplified amplitudes IBl and IB2 



(PHOTOS-type and Ginsberg- type) . We hope that our 
results can help construct more reliable MCs for the sim- 
ulation of B irevj decays, or at least address the sys- 
tematic uncertainty to be assigned to the present MC re- 
sults. The codes used to generate the distributions given 
in this section are available from the authors upon re- 
quest. 



VI. CONCLUSIONS 

We derived in this paper a new factorization relation 
for the exclusive semileptonic radiative B — > jniiy de- 
cays, valid in the kinematical region with a soft pion and 
an energetic photon. In this region, the form factors sat- 
isfy the factorization relation written in schematic form 
in Eq. The factorization relation is easily general- 
ized to any process B ^Miv, with M = p,uj, - ■ ■ a 
light soft meson. 

In all cases, this relation contains as input parameters 
the nonperturbative matrix elements of two soft opera- 
tors for the B M transition. These matrix elements 
are the B physics analog of the generalized parton distri- 
bution functions encountered in the factorization relation 
for the off-diagonal 7*A^ — > N'j transition, or 'deeply 
virtual Compton scattering' (DVCS) 

For sufficiently slow pions, the _B — > tt soft function 
can be computed using chiral perturbation theory meth- 
ods. Previous applications of chiral symmetry to con- 
strain GPDs for light quark systems have been given in 
|3^ . At leading order in A/mt and in the chiral expan- 
sion, the only nonperturbative quantity required is the B 
meson light-cone wave function. 

We have presented phenomenological applications of 
our factorization formula, with two main goals in mind. 
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On one hand, we have given factorization predictions for 
radiative observables (such as photon spectrum and in- 
tegrated rate) that can be tested as increasing statistics 
becomes available at the B factories. Such a test would 
provide an important validation of the theoretical tools 
used here. On the other hand, we have compared our 
results with simplified approaches to radiative hadronic 
decays, relying on universal properties of radiative ampli- 
tudes following from gauge invariance, not taking into ac- 
count hadronic strucure effects. These simplified ampli- 
tudes are currently implemented in the MC simulations 
used in experimental analysis. In the kinematical region 
accessible to our theoretical tools (energetic photon and 
soft pion in the B rest frame), we find a large discrepancy 
between our factorization results and PHOTOS-type or 
Ginsberg- type results. We hope our work will provide a 
theoretical basis to construct improved MCs for the sim- 
ulation of B — > ireDj, which is a crucial ingredient in 
a high precision measurement of the semileptonic mode 
B — > TTeD. 
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APPENDIX A: LOW'S AMPLITUDE 



In this Appendix we report the result for A{B{p) 
'K{p')e~{pe)v{pi,)^{q,e)) predicted at small q by Low's 
theorem 8] . This serves as basis of simplified treatments 
of the radiative process 0, , and is used in the text 
for a comparison with the results of our paper following 
from the factorization relation. 

Assuming E-^ ^ A = m^* — friB, one finds for the 
neutral mode decay: 



A 



ubuipe) 



e* -p' , e* ■ pe 



eGp 

\/2 ' LV q-p' ' q-pe ' 2(7-p£ 



V2 



■Vubu{pe) 2p-D{p'){2i)- 



df+ dh\ 



dt^J 



T"{p,p')\v{p.) 
(1 - 75) v{p^) , 



(Al) 



where the form factors f±{t'^) appear in the parameterization of the matrix element F^{p,p') given in Eq. ((HJ and are 
defined as 



F^,{p,p') = f+{t^){p + p')^ + f-{t'){p-p')^ 
with t^p-p'. We define also /i,2(t^) = /+(t^) ± f-it^), and 

T'{p,p') = m+{e}~mj,{e}]{i-j,) 

D{v) = ^—^q^£* . 
q ■ V 



(A2) 

(A3) 
(A4) 



r 



For reference, we quote here also the definition of the 
B ^ TT tensor form factor, used in Sec. IV 

{7T{p')\qta^,b\B{p)) = fT{t^){p^pl - Pup') (A5) 



In the limit ni^ — » Low's amplitude leads to the 
following expressions for the invariant form factors ap- 
pearing in Eq. jSjl: 



Ui 
^2,3,4 
^1,2,3,4 



Q-P 



l-2q-W 



1 d.f+- 







(A6) 
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Similarly, for the charged mode decay we find: 

3 eGp 



A 



Low 



V2 
eG 



■VubU{pe) 



s-p , e -Pe 



q-p 



q -Pe 'iq-Pe 



T {p,p')^v{p^) 

''^Vub u(pe) [2p' • D{p) (2^'^ + m,|^) (1 ~ 75)] v{p.) , 



(A7) 



where 

T-{p,p') - [2fU{t') + mj,{t^)]{l - 75) . (A8) 
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